Direct Search for Low Mass Dark Matter Particles with CCDs 
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A direct dark matter search is performed using fuUy-depleted high-resistivity CCD detectors. 
Due to their low electronic readout noise (R.M.S. ^7 eV) these devices operate with a very low 
detection threshold of 40 eV, making the search for dark matter particles with low masses 5 
GeV) possible. The results of an engineering run performed in a shallow underground site are 
presented, demonstrating the potential of this technology in the low mass region. 



PACS numbers: 93.35. +d, 95.55.Aq 
I. INTRODUCTION 

There have been several direct-detection experiments 
searching for dark matter (DM) performed in recent 
years, and several more in development [1 . Most of these 
experiments have been optimized for detecting the elas- 
tic scattering of DM particles with masses larger than 
50 GeV, concentrating on the most natural region of the 
Weakly Interactive Massive Particle (WIMP) parameter 
space such that DM corresponds to the lightest super- 
symmetric particle 0. Detection thresholds of a few 
keV are typical for such high mass DM searches. Re- 
cent results from experiments that operate with a lower 
threshold [3, 4 have presented hints of a DM signal at 
low energies, the most notable being DAM A/Libra which 
claims a high-significance detection of an annual modu- 
lation ^ . More recently, the CoGeNT Collaboration has 
seen a possible hint for a low mass DM signal using a low 
threshold Ge detector [4 . Motivated by these results, 
new theoretical interpretations have been developed for 
which low thresholds are needed to directly detect DM 
[5HTTj . In these models the DM particles are either much 
lighter than 50 GeV or present a non-elastic interaction 
with the detector nuclei [12]. Because of noise in detec- 
tor electronics and efficiency issues, most experiments are 
not capable of significantly lowering the detection thresh- 
old to probe the new models. For this reason there is an 
effort by the experimental community to develop tech- 
niques that could detect low energy signals from DM 
[4] [13]. This work presents a DM search using Charge- 
Coupled Devices (CCDs) using a detection threshold of 
40 eV electron equivalent energy (eVee). 

CCDs used for astronomical imaging and spectroscopy 
commonly achieve readout noise levels of 2 e~ R.M.S. , 
equivalent to 7.2 eV of ionizing energy in Silicon. Though 
this allows for a very low detection threshold, these detec- 
tors have not been considered for DM searches because 



of their very low fiducial mass. The recent development 
of thick, fully-depleted CCDs, ten times more massive 
than conventional CCDs, has changed the situation. The 
Dark Matter in CCDs (DAMIC) experiment is the first 
DM search to exploit this technology. 



II. HIGH RESISTIVITY CCD DETECTORS 

In an effort to increase the near-IR photon detection ef- 
ficiency of CCDs, a group at Lawrence Berkeley National 
Laboratory (LBNL) developed detectors with a depletion 
region up to 300 jiia thick [14 . These devices are fab- 
ricated with high resistivity silicon (~ lOkl^ cm). Due 
to their higher efficiencies for detecting photons with red 
and near infrared wavelengths, these devices have been 
selected for astronomical instrumentation [15] [16], in par- 
ticular, the Dark Energy Camera (DECam) [I7l[l8]. The 
DAMIC experiment uses DECam CCDs to search for 
DM. 

The DECam CCDs [1^ IMSIl are back- 

illuminated, p-channel CCDs thinned to 250 /im and bi- 
ased with 40 V from the back side to achieve full deple- 
tion. The positively-charged holes produced in the deple- 
tion region are stored in the buried channels, a few /im 
away from the the gate electrodes. Charge produced near 
the back surface must travel the full thickness of the de- 
vice to reach the potential well. During this transit inside 
the depletion region, a hole could also move in the direc- 
tion perpendicular to the pixel boundaries. This effect, 
called charge diffusion, must be regulated to avoid a sig- 
nificant degradation in image quality. The CCDs used in 
most astronomical instruments are thinned to <40 /im to 
reduce charge diffusion. For the DECam CCDs, a sub- 
strate voltage is applied to the back surface to control 
diffusion and obtain acceptable image quality. 

The two main features that make the DECam CCDs 
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FIG. 1: Four muons tracks are clearly reconstructed on this 
image. There is a significant width different between the two 
ends of the track. The thicker end corresponds to the region of 
the track in the back of the CCD, with maximum diffusion. 
The thinner end corresponds to the front of the CCD. The 
smaller circular hits are diffusion-limited hits as expected for 
nuclear recoils. The line is 100 pixels long and was added to 
give a sense of scale. 



good candidates for a direct dark matter search are their 
thickness, which allows the CCD to have significant mass, 
and the low electronic readout noise, which permits a 
very low energy threshold for the ionization signal pro- 
duced by nuclear recoils. 

DAMIC uses a single rectangular CCD readout with a 
Monsoon controller [22 . The detectors have 4.2 million 
pixels with dimensions 15 /im x 15 /am (active mass of 
0.5 g), and are read by two amplifiers in parallel. The 
detectors have an output stage with an electronic gain of 
~ 2.5 /iV/e. The signal is digitized after correlated dou- 
ble sampling (CDS) and the noise performance depends 
on the readout speed. The best noise, a < 2e~ (R.M.S.) 
per pixel, was achieved with readout times of 50 /is per 
pixel [19H2I1E3]. 



III. CCD CALIBRATION: X-RAYS AND 
NUCLEAR RECOILS 

X-rays are commonly used in the energy calibration 
of CCD detectors [24 . X-rays from ^^Fe penetrate only 
~20 jam into the silicon before producing charge pairs ac- 
cording to the well known conversion factor of 3.64 eV/e- 
[24 . As a result of this process X-rays produce hits in 
the detector with size limited by charge diffusion. In a 
back illuminated DECam CCD the size of the 5.9 keV X- 
rays is ~ 7.5 /im R.M.S. [25]. Using processing tools de- 
veloped for astronomical aplications [26 , the X-ray hits 
are identified in the CCD images as reconstructed charge 
clusters with a size limited by diffusion. Fig. [l] shows 
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FIG. 2: Spectrum obtained for the reconstructed X-ray hits 
in an ^^Fe exposure of a DECam CCD. The arrows mark the 
direct X-rays from the source: Ka=5.9 keV and K/3=6.5 keV, 
their escape lines at 4.2 and 4.8 keV 24 , and the Si X-ray at 
1.7 keV. The factor 3.64 eV/e is used to convert from charge 
to ionization energy. The feature at 7.6 keV is consistent with 
pixels that hit by a 5.9 keV and a 1.7 keV X-ray on the same 
exposure, consistent w 



two examples of diffusion-limited hits. The energy of the 
X-ray hit is proportional to the charge collected in all 
the pixels that belong to this cluster. An example of the 
energy spectrum measured for an ^^Fe X-ray exposure in 
a DECam CCD is shown in Fig. [2] 

Figure [2] provides a good calibration for X-rays in sili- 
con. Since the ionization efficiency for nuclear recoils dif- 
fers from the ionization efficiency for X-rays, the results 
are not directly applicable to nuclear recoils expected 
from DM interactions. The ratio between the ioniza- 
tion efficiency for nuclear recoils and ionization efficiency 
for electron recoils is usually referred to as the quench- 
ing factor (Q). The quenching factor has been measured 
in Si for recoil energies above 4 keV [27], showing good 
agreement with the Lindhard model [28l, ^29^ . For recoil 
energies less than 4 keV, the quenching factor becomes 
increasingly energy-dependent and no previous measure- 
ments are available, as shown in FigjS] In this work, the 
Lindhard model [28l is used for converting visible (or 
electron equivalent) energy to recoil energy. We test this 
assumption by comparison to neutron data with a known 
energy spectrum. 

In order to study Q, the DAMIC detector was exposed 
to a ^^^Cf neutron source (see Fig|4|. Nuclear recoils pro- 
duced by neutrons, similar to those expected from DM in- 
teractions, give diffusion-limited hits analogous to those 
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FIG. 3: Existing quenching factor measurements in Si com- 
pared with the Linhard theory (soUd hne). No previous data 
exist for recoil energies below 4 keV ^28l [29] . 



of X-rays. Thus, nuclear recoils in CCD images are iden- 
tified by selecting diffusion-limited hits. The expected 
response of the detector was calculated in three steps. 
First neutrons were generated according to the known 
spectrum of the ^^^Cf source and GEANT4[30 was used 
to simulated the passage of these neutrons through the 
vacuum vessel wall housing the CCD detectors ( 1cm Al 
wall). The resulting neutron energy spectrum was then 
used to calculate the expected nuclear recoil in the Silicon 
CCD using Ref.[27^. Finally, using the ionization yield 
from the Lindhard theory, the expected observable en- 
ergy on the CCD was calculated. The results of the sim- 
ulation compared with the spectrum for diffusion-limited 
hits in data are shown in Fig. [5] The data shows a bump 
at ^1.7 keVee consistent with Si excitations, but deviates 
from Lindhard model at energies below 1.5 keV. While 
the data appears to indicate a weaker energy dependence 
of the quenching factor below 1.5 keV, we use the Lind- 
hard model in order to produce more conservative limits. 
For energies lower than 0.5 keV the detection efficiency 
for nuclear recoils has a strong energy dependence (see 
Figj9| and the comparison is no longer valid. This issue 
will be investigated in future work. 

In Fig js] the behavior of the ionization efficiency of nu- 
clear recoils is degenerate with the selection efficiency. 
The spectrum is also contaminated with electron recoils 
produced from the gammas generated in the ^^^Cf souce. 
The comparison between data and simulation is also in- 
fluenced by the neutron input spectrum and the simu- 
lation geometry. For these two reasons the analysis dis- 



FIG. 4: Image resulting from an exposure of a DEGam GGD 
to a ^^^Gf neutron source. The total width of the image 
corresponds to 1000 pixels. The smaller dots represent the 
diffusion- limited hits, the trails correspond to scattered elec- 
trons and there is one bigger circular cluster of charge that 
corresponds to an alpha particle. 
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FIG. 5: Reconstructed electron equivalent energy spectrum 
for ^^^Gf exposures. The data is consistent with expectations 
from Lindhard theory (red). The expectations for an energy 
independent quenching factor are also shown for comparison 
(green). 

cussed above, and summarized in FigjS) does not consti- 
tute a measurement of Q. It should be interpreted as 
a comparison between the data and Lindhard model as- 
suming a constant detection efficiency. It is presented 
here to illustrate the effect of energy dependent Q in the 
recoil spectrum and to motivate the need for a full cali- 
bration of nuclear recoils in Si at low energies. 



IV. SELECTION OF DARK MATTER 
CANDIDATE EVENTS 

Three selection cuts are used to separate the dark mat- 
ter candidates in our images from background and noise 
hits. The first step requires the total energy deposited 
to be larger than 0.04 keVee. This cut is mostly used to 
suppress the noise on the readout of the CCD detector. 
In order to select diffusion-limited hits produced in the 
bulk of the CCD and not near the front or back surfaces, 
we impose additional two selection described below. 
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X-rays produce diffusion-limited hits similar to nuclear 
recoils in silicon, becoming an important background for 
a dark matter search. However, unlike nuclear recoils, the 
low energy X-ray photons only penetrate a few microns 
into the detector, producing charge very close to the back 
or front surfaces of the CCD. When performing a DM 
search it is convenient to reject hits produced near the 
front and back surfaces of the CCD. 

Hits on the front surface will have negligible diffusion 
because they are produced next to the gates. Hits on the 
back surface have the maximum diffusion since the path 
is the longest to the electrodes. The dependence of diffu- 
sion with depth can be clearly seen in tracks identified as 
cosmic ray muons in Fig. [l] The diffusion measured for 
each hit becomes indicative of the depth of the hit inside 
the detector. 

For each hit reconstructed on the CCD image the first 
order moments of the charge distribution along each axis 
are calculated, namely x and y. These values are an 
estimation of the real position of the hit on the CCD 
coordinates, x and y. The second order moments along 
both axes are also calculated as and a^, providing 
an estimation for the charge spread and a^. These 
four estimators are the basis of the selection cuts used 
to eliminate the events near the front and back surface. 
In the following discussion units of pixel and pixel^ are 
implicit when presenting values for x and a^. 

Unfortunately, because of the pixelation of the detec- 
tors there is a non trivial relation between the real val- 
ues of X, y, cr^, ay and the estimators ^, a^, discussed 
above. For example, events close to the boundary of the 
pixel necessarily have > a^. At the same time, events 
with small charge spread tend to have x biased to- 
wards the center of a pixel. Also due to pixelation, there 
is a strong correlation between and given by 

al > 0.25 - x^ for 0.0 < f < 0.5 and 

(1) 

al > 0.25 - (x - 1)^ for 0.5 < x > 1.0 . 

With X = 0.0 corresponding to the edge of the pixel and 
X = 0.5 the center of the pixel. A similar relation applies 
for the y-ax.is. 

A simulation was performed to illustrate the pixelation 
effects on x and . The simulation consisted of gener- 
ating point-like charge distributions in a 250 /im thick 
CCD with 15x15 jam^ pixels. Lateral diffusion for each 
hit was calculated as a function of the distance of the hit 
to the front surface (maximum diffusion was 7.5 /im for 
hits 250 /im away from the front). The results shown in 
Fig. [g] clearly indicate that the larger is measured for 
events on the back of the detector, and that events on 
the front surface are very close to the curve defined by 
Eq.g. 

The results for x and with X-rays data on the front 
and back of the detector and shown in FigjTj together 
with with data from a ^^^Cf neutron source. The data 
shows consistency with the simulation. 

Based on the discussion above, the second step for the 
selection of DM candidate hits consists of requiring the 




FIG. 6: Results of a simulation of charge diffusion inside a 
CCD. The horizontal axis shows the position of the hit inside 
the pixel, with x=0.5 corresponding to the center of the pixel. 
The vertical axis shows the size of the reconstructed hit, 
in units of pixeP as defined in the text. The red points cor- 
respond to hits generated in the first 50 /im adjacent to the 
front of the detector, the green points correspond to the hits 
generated less than 10 /im away to the back of the CCD. The 
black points correspond to hits generated between 50 and 200 
/im from the front of the detector. The simulation shows that 
the selection criteria discussed in Eq.([2]) and Eq.(|3| rejects 
most events next to the front and back surfaces of the CCD. 

charge to be distributed in more than one pixel. This is 
used to reject X-rays in the front surface of the detector. 

The third selection cut is designed to keep hits in the 
bulk of the CCD while rejecting events close to the front 
and back surfaces of the detector. This is done in two 
steps, the requirement 

al - (0.25 - f ^)) >0.05 for 0.0 < f < 0.5 and 
al - (0.25 - (1 - xf) >0.05 for 0.5 < x < 1.0, 

is used to ensure that the selected events are away from 
the front of the detector. To select events away from the 
back of the detector one additional condition is added 

al < 0.28 (3) 

Identical cuts are applied to the ?/-axis variables. 

The efficiency for nuclear recoils passing the selection 
is estimated using ^^^Cf exposures. The reconstructed 
spectrum with and without the selection is shown in 
Figjs] The raw spectrum (without size selection) is con- 
taminated by low energy photons and electrons that pro- 
duce a steep raise below 0.5 keV. This contamination is 
also evident in the image shown in Figj4j Nuclear recoils 
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FIG. 7: Results of x and al for 5.9 keV X -rays on front (red) 
and back (green) of a CCD. The black data corresponds to 
hits from the same ^^^Cf source shown in Fig|4] As expected 
from simulations the back illuminated X-rays produce hits 
with large a^, while the front illuminated X-rays are very 
close to the boundary defined by Eq^ 
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FIG. 8: Measured spectrum for ^^^Cf . The upper line shown 
the raw data without applying any selection criteria, the lower 
hatched histogram shows the result after selecting hits with 
intermediate size. The raw data is fitted to a linear function 
above 0.5 keV. 



from neutrons emitted by ^^^Cf are expected to produce 
an approximately flat spectrum at low energies. We de- 
termine the efficiency of the selection by taking the ratio 
of the selected hits to a linear fit to the data before selec- 
tion, in the region 0.5 keV to 10 keV where the spectrum 
is the most consistent with a flat line (Fig. |8|. The result 
is shown in Figj9] 
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FIG. 9: Efficiency for selecting nuclear recoils with interme- 
diate size as calculated using the data in Fig|8] 




FIG. 10: Schematic of DAMIC vacuum vessel and lead shield. 



V. TEST IN SHALLOW UNDERGROUND SITE 

We performed the first underground test for the use of 
CCDs in a direct DM search with a single 0.5 g CCD as 
described above. A low background CCD package was 
designed for this test, consisting of a readout board built 
on aluminum nitrade (AIN) wire-bonded to the CCD 
with some additional AIN spacers for mechanical sup- 
port. AIN is typically used in CCD packaging due to 
its mechanical and thermal properties and was measured 
to have low concentrations of radio-isotopes. The CCD 
package was installed inside a copper box and cooled 
down to -150C to reduce dark current intrinsic to the 
CCD. The copper box was also cooled to eliminate the 
infrared radiation impinging on the detector from warm 
surfaces. A closed cycle helium gas refrigerator [31] was 
used to maintain the low temperature. The detector was 
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TABLE L Number of events passing selection cuts for the 107 
g-day of data shown in FigjlT] two energy bins are shown: 
(1) from 0.04 keV to 5 keV , (2) from 0.04 keV to 15 keV. 



cut 


bin (1) 


bin (2) 


1) E > 0.04 keV 


81754 


102469 


2) npixel > 1 


26971 


45353 


3) hit size 


433 


5529 



connected with a passive readout cable and the active 
electronic components (preamplifiers) were located out- 
side a 6" lead shield. The detector package was housed in 
a cylindrical vacuum vessel fabricated with oxygen-free 
copper, and maintained at 10~^ Torr by continuously 
running a turbo molecular pump as in FigfTo) The ap- 
paratus was installed 350' underground in the NuMI [32] 
near-detector hall at Fermilab, where other DM exper- 
iments have already been performed [33 . The system 
operates without the need of any human intervention, 
and was programmed to collect exposures every 40,000 
seconds. The experiment operated for 11 months start- 
ing June 2010, and accumulating a total exposure of 107 
g-days. The observed spectrum in shown in Fig. 11 The 
spectrum shows several X-ray peaks coming from known 
components the apparatus. It also shows three peaks at 
11 keV, 12 keV and 14 keV for which the source has not 
been identified. The isotopic composition of the commer- 
cial cable and high density connector used for the CCD 
have not been fully characterized, and are the most likely 
source of these unidentified features. These X-ray peaks 
limit the DAMIC reach in a DM search, and point to the 
need of improving the isotopic control of the apparatus. 

The low energy nuclear recoil candidates were selected 
as described above and the resulting spectrum is shown 
Given the shallow depth of the underground 
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in Fig. 

site, we expect that most of the events in the spectrum 
correspond to neutrons produced from cosmic ray muons 
hitting the rock or the lead shield. Because of the lack of 
timing information from CCDs due to the long exposure, 
DAMIC cannot benefit from an active muon veto as is 
commonly used in other DM searches. The number of 
events passing each selection criteria are shown in Table 

ID 



VI. RESULTS AND CONCLUSION 

Standard techniques described in Ref. [27^ were used 
to interpret these results as a cross section limit for spin- 
independent DM interactions. We assume a local WIMP 
density of 0.3 GeV/cm, dispersion velocity for the halo 
of 230 km/s, earth velocity of 244 km/s and escape ve- 
locity of 650 km/sec. The Lindhard model was used to 
obtain recoil energies as discussed above. Other possible 
combinations of parameters could have been used, in par- 
ticular, the choice of escape velocity might be considered 
optimistic for this analysis. However, these choices were 



made to compare directly with recent publications from 
experiments looking for low mass dark matter particles 
[4. 

The optimal interval method j34ll37] was used for de- 
termining the upper limit on DM cross section as a func- 
tion of mass. The energy range used for this analysis 



is between 0.05 keVee and 2 keVee as shown in Fig 12 
The resulting 90% C.L. limits are shown in Fig. 13 
and constitute the new best limit for dark matter par- 
ticles of masses below 4 GeV. This corresponds to an im- 
provement over the limits produced by CRESST-I [38], 
which were achieved using sapphire cryogenic detectors 
with a threshold of 600 eV and exposure of 1.5 kg-day. 
The region consistent with a DM interpretation of the 
DAM A/LIBRA signal is shown as a shaded area in Fig. 
[T3| , along with recent results from the CoGeNT collabo- 
ration interpreted as a DM signal [4l[9]. These two results 
make the low mass region very interesting and clearly in 
need of further studies. 

The test at the shallow site presented here shows an 
average background rate of 600 cpd/keV (see Fig. 11). 
This initial setup did not have a shield for neutrons com- 
ing from the rock. Because of the lack of timing infor- 
mation on CCDs, it is not possible to take advantage of 
an active muon veto to remove cosmogenic neutrons. For 
these reasons the largest contribution to the background 
comes from neutrons. The next step for this technology 
becomes then the development of a CCD experiment in a 
deeper site (where cosmogenic neutrons become negligi- 
ble) and with an efficient neutron shield. A CCD search 
for DM with a background of 10 cpd/keV and a thresh- 
old of 40 eVee would improve the limit presented here 
by a factor of 60, and could provide a new probe for 
the low mass region where hints of DM have been ob- 
served by other experiments. The DAMIC collaboration 
is currently considering an installation of the experiment 
a SNOLAB to reach this goal. 

Recently we have demonstrated that 0.2 e- R.M.S. 
readout noise can be achieved with a new type of 
CCD detector fabricated by LBNL using the same high- 
resistivity Si technology [42 , but with a non-destructive 
output circuit enabling multiple readouts to reduce noise, 
a so-called "skipper CCD" [24]. This corresponds to a 
readout noise 0.72 eV and a reduction of a factor of ten 
compared to the results presented here. A DM search 
with the skipper CCD would improve significantly the 
reach at low energies and is being considered as part of 
a DAMIC upgrade. 
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FIG. 11: Energy spectrum for nuclear recoil candidates mea- 
sured in a 107 g-day exposure for DAMIC. The upper his- 
togram shows the spectrum for all hits, the solid line just 
below is the result of eliminating the single pixel hits, and 
the hatched histogram corresponds to the events passing the 
selection cuts described in the text. The 8.0 keV peak can be 
attributed to Cu X-rays from the material surrounding 
the CCD. The 5.9 keV and 6.4 keV peaks can be attributed 
to the Mn Kq, X-rays most likely from cosmogenic ^^Fe. The 
9.6 keV peak matches the Au La line. The peaks at 11 keV, 
12 keV and 14.0 keV match Ge, Br and Sr X-rays, but these 
elements can not be attributed to any known component in 
the DAMIC apparatus. 
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FIG. 12: Low energy spectrum of events passing the selection 
cuts. The vertical lines at 0.04 keV and 2 keV show the energy 
range used for the DM analysis. 
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